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Abstract. We model the immune surveillance of a pathogen which passes through n immunologically distinct 
stages. The biological parameters of this system induce a partial order on the stages, and this, in turn, determines 

f^^ ' which stages will be subject to immune regulation. This corresponds to the system's unique asymptotically stable 

^vj ' fixed point. 

1. Introduction 

" ' ' Pathogens that traverse different stages during their Hfe cycle or during an infection process have been 

studied since the late nineteenth century. The most prominent genus is Plasmodium, causer of Malaria, 
Py ' a pathology first described by Hippocrates in the fourth century B.C., which is still a major public health 

Ph ■ problem in many continents yt]. Other important examples are Trypanosoma 12511 and Epstein-Barr virus 

Q ' (EBV) which is a member of the family of herpes viruses. Indeed, it seems likely that this entire family 

•^ . uses multiple stages 12311 . Our focus is on Epstein-Barr virus, which is known to cycle through at least 



o 



S 



four different stages during infection within the human body 112211 . (For a quick sketch of EBV biology, see 
O^. Section|9]|. One remarkable characteristic of infections with many of such pathogens is life-long persistent 

infection ^, il, HI, H, Q. 



Our main goal in this article is to study the properties of the immune response to such a pathogen us- 
^ ' ing mathematical modeling. Mathematical approaches to studying host-pathogen interactions have steadily 



Cn , increased in the last 4 decades. For entry into the corresponding body of literature, we recommend 11611 

iL: ' il8ll and 12611 . What to include in such a mathematical model depends not only on the host and pathogen in 

^— N . question, but also the taste and purpose of the modeler. Thus models vary in scope and complexity. A simple 

model might include uninfected tissue, infected tissue and an immune response. More complex models can 
involve such factors as the details of T-cell activation, competing immune responses to multiple epitopes and 
p^ ' the mutation, selection and consequent evolution of the pathogen. Much of this literature is driven by the 

urgent need to understand the host-pathogen dynamics of HIV infection. (A Pubmed search on the terms 
"HIV" and "mathematical model" produces close to 400 hits.) Our model is motivated by the broad outlines 
of EBV biology which it both simplifies and generalizes. 
/^ ' In the mathematical model we propose herein, the pathogen traverses a cycle of immunologically distinct 

stages during infection of the host. On the pathogen side, our model includes the natural proliferation or 
decay of the pathogen population at each stage, the rate at which each stage is lost to produce (or become) the 
next stage and any gain factor involved in this process. (For example, one lytically infected cell may produce 
^ 10^ free virions.) On the host side, it includes a distinct immune response to each stage, proliferation of 
this response due to encounter with the pathogen, the ability of this response to kill or disable the pathogen, 
and the decay of the response in the absence of the pathogen. Omitted from this model is the supply of 
uninfected tissue which supports the cycle of infected stages. This omission is only likely to be important 
when this supply is a limiting factor This could happen when the immune response is weak or non-existent. 
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as when a naive host is first exposed or in the case of an immunosuppressed host. As a consequence, this 
model may be more realistic in describing the long-term behavior of persistent infection when this is not 
tissue-limited. 

We will see that when the host and pathogen stand in this relation to each other, there is a unique stable 
accommodation between them. This appears in the model as the unique biologically meaningful asymptot- 
ically stable fixed point. In particular, it determines exactly which stages will fall under immune regulation 
and determines the levels of pathogen population and host response at which this mutual accommodation 
occurs. 

In Section|2]we will present the model's equations and the range of parameter values under consideration. 
We will start by considering parameter values such that the pathogen must traverse all stages in order to 
establish infection. In particular, no stage is capable of proliferating and establishing infection independent 
of the others. In addition, we will restrict our attention to generic values of the parameters and will exclude 
certain parameter sets of measure zero. 

The central topic of this article will be the analysis of the steady state behavior of our model. To this end, 
we start by calculating necessary and sufficient conditions for the pathogen to establish infection. We view 
this both in terms of the basic reproductive constant, Rq, and in terms of the linear stability analysis of the 
uninfected equilibrium at which all stages of the pathogen as well as their corresponding immune responses 
become extinct. 

We continue the steady state analysis by determining further fixed points of the model's equations. As 
we will see, for generic parameter values, this system has exactly 2" fixed points, (where n is the number of 
stages), depending on which stages are regulated by the immune response. In general, not all of these are 
biologically relevant since a given pattern of regulated and unregulated stages can give negative population 
levels at steady state. We give a criterion for detecting the biologically meaningful ones. 

If the pathogen is able to establish infection, [Rq > 1), the parameters induce a partial order on the 
pathogen's stages. We say stage j starves stage k (in symbol j >- k) if immune regulation at stage j deprives 
stage k of sufficient population to support immune regulation. A stage k is called starvable if there is another 
stage j such that j >~ k. If no such j exists, k is called imstarvable. One of our main results is the fact 
that, generically, the system has a unique asymptotically stable fixed point, namely, the one at which all 
unstarvable stages are regulated and all starvable stages are unregulated. 

The portrait we have just drawn holds for generic values of the parameter set. It turns out that there 
are highly non-generic parameter sets for which the system becomes a Lotka- Volterra predator-prey system. 
While mathematically interesting, this is unlikely to occur in nature. The condition that the parameters must 
obey has co-dimension « — 1 . In particular, such parameters form a set of measure zero. 

Finally, we show how to generalize the previous results to the case where the pathogen has individual 
stages which are capable of establishing infection independently of the remaining stages. 

This article is organized as follows; In Section|2] we introduce the model and discuss its parameters. In 
Section[3] we examine the basic properties of the system's fixed points. In Sections HI |5] and |6] we embark 
on the detailed analytical study of the model's equilibrium behavior. One of the main tools used is the 
partial order on stages, defined in Section|5] In Section]?] we present the highly non-generic Lotka- Volterra 
scenario. Finally, in Section [8] we consider self-establishing stages. We close in Section |9] with discussing 
some of the model assumptions and consequences for the characteristics of the immune response to a multi- 
stage pathogen. To provide a more concrete contrast between the model and reality, we focus on EBV 
infection biology. 
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2. Formulation of the model 



We model the interaction between a host and a pathogen which traverses « e N, « > 1 different stages 
during its life-cycle. We treat these stages as immunologically distinct, that is, the host produces a separate 
immune response to each of them. Our model is given by 

(*) Sj = Fj{S, T) = rj^ifj^iSj-i - ajSj - fjSj - pjSjTj 

tj=Gj{S,T)=cjSjTj-bTj. 

The subscripts are taken modulo n. We use Sj and 7} for the pathogen population and immune responsqj 
at each stage. These are assumed non-negative. The parameters represent the following processes: 

• Oj is the decay rate of stage Sj. If Oj is negative, this state proliferates. 

• fj is the rate at which stage Sj is lost to become (or produce) stage Sj+\. 

• rj is an amplification factor in the process by which stage Sj becomes (or produces) stage Sj+i. For 
example, the loss of one lytically infected cell may produce rj = 10^ free virus. 

• pj represents the efficacy of the immune response 7} in killing infected stage Sj. 

• Cj is the antigenicity of stage Sj, i.e., its efficacy in inducing proliferation of immune response Tj. 
9 b is the natural death rate of the response Tj. We assume it is the same for all stages. 

We will refer collectively to a choice of values for these parameters as n. We will assume that with the 
possible exception ofay, all of these are positive. In Sections[3]through|7] we will assume that Oj + fj > 0. 
Thus, while some stages may proliferate, none does this as fast as it is consumed in producing the next 
stage. In Section [8] we will lift this assumption and generalize our results to the case in which there are 
self-establishing stages. We will also assume that the values of these parameters are generic. Accordingly, 
at several places we will disregard certain sets of parameters which have measure zero. (See Equations (O 
and ©.) 

The model's equations can be simplified by using the following change of coordinates 



Tj = PjTj 



which produces the equations 



Sj = Fj{S,T) = rj^ifj^iSj-i-ajSj-fjSj-SjTj 
(**) tj = Gj(S,T) = SjTj - bTj 



rj = rj 

a 



We will henceforth assume our equations are in this form and omit the bars. When we refer to a fixed point, 
we will mean a fixed point of system (Q. The notations {Sq,...,SI_i,Tq,...,T*_^) or {S*,T*) will refer to 
such a fixed point. We will say that a point is biologically meaningful if its components are non-negative. 

We will adopt the following notational conventions. Sets such [j,k] and [j,k) are to be taken cyclically. 
That is to say, if j < k, then [j, k] — {j,... , k}, while if j > k, [j,k] = {j,... ,n — l,0,...,k}. We take [j,j) 
to be the empty product so that [j,j) = 1 . We abuse notation by taking [0, n) = {0, . . . , n — 1 }. 



The notation Tj is motivated by the T-cell response. However, it may equally well refer to humoral response. 
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3. First properties of the fixed points 
The fixed points are defined by the system of polynomial equations 

(1) Sj = Fj{S*,T*)^rj^ifj^iS*_,-S*{aj+fj + T*)=0 

(2) tj = GjiS*,T*)^{S*-b)T* =0 

where j = Q,...,n — l and indices outside the interval [0,n — 1] are understood modulo n. 

In this section, we give formulas that determine the population values at fixed points of system (F*T l. We 
will see that if stage k is unregulated, S*j is determined as a "follow-on" population from the previous stage. 
This gives rise to "follow-on factors" (defined below) which will play a key role in this paper. This, in turn, 
will lead us to the notion of generic parameter values, and we will see that for generic parameter values, 
system i**[ has exactly 2" fixed points. 

This section is organized as a series of observations. These follow from easy computations which we 
omit when they would clutter rather than clarify. 

We start by defining these follow-on factors 

Mj := 'J^ 

Mjk-^ n ^( 

(W-.k) 

Ro := n Mj- 

je[0,n) 

which we will use throughout this article. We will take Mjj to be the empty product. Given a fixed point 
{S*,T*), if r* ^ 0, we will say that j is regulated. Otherwise j is unregulated. We define 

Reg(5*,r) = {;|7;*^0} 
Unreg(5*,r) = O|7;*=0} 

1) If there is a j such that S* = then {S* ,T*) = (0,0). By © if 5* = 0, then S*_^ =0. By © if 
S* == 0, then T* = 0. Continuing in this way S* ^T* =0 for all ; G [0,n). 

2) Ifj&Rcg{S*,T*),thenS*^b. 

3) //;■ G Unreg(5* ,T*),thenT*= 0. 

4) //;■ + 1 G Unreg(5*.r*), then S*^^ = S*Mj. This follows quickly from ©. 

5) ^[y + l;^] cUnreg{S*,T*), then SI ^ S*Miit. This follows by induction on the previous observa- 
tion. 

6) IfjeRegiS*,T*)and[j + l,k] cVmeg{S*,T*), then S; =bMji,. 

7) IfRo ^ 1, and {S*,T*) ^ (0,0), then Reg{S*,T*) ^ 0. To see this, notice that by our first observa- 
tion, we must have 5* y^ for all ;■ G [0 , « ) . But if Reg (5* , r * ) = 0, we must have 5o = /?o5o . 

Before proceeding, two points are worth considering. The first is that we have just introduced our first 
condition for a generic parameter set n, that is, 

(3) Ro= n ^J^^ 

ie[o.«) 

The second is that the S population at an unregulated stage depends on follow-on constants and the previous 
regulated stage. Observation 7 guarantees that at a non-trivial fixed point, there always is such a stage. Ac- 
cordingly, given k, we will define h/^ to be this stage, that is, h/^ is the unique stage such that /i^ G Reg(5'*, T*) 
and [hii + l,k) C Unreg(5'*,r*). As we will see, this notation is also useful in the case k G Reg(5'*,r*). 
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8) If j e Regis*, T*), then T* = '"'~'/^'' ^}_i - (aj + fj). This follows from ©■ 

9) Ifj e Regis*, T*), then T* = ry_i/y_iM/, .^^i - (fl/ +/,). This follows from the fact that 5*_i = 
SlMh.j^i and Sl^=b. 

10) T* > if and only ifMh^j > 1. 
This latter will be important in Section|5]and leads to our second genericity requirement. We take 

(4) M,-^^lfor;V^- 

This may seem overly broad since we have only used Mji^ ^ 1 for fk G Reg(5'*, T*), with j = hj^. However, 
we wish to avoid making the genericity of k depend on the fixed point in question. The sets in question ai-e 
all of co-dimension 1 and hence have measure 0. 

We have now arrived at the place where all of the populations are determined by the parameters and the 
pattern of regulation. We record these as follows: 

(5) S* = b and T* = r,_i/,_iM,,^,_i - (a,- +/,) V ; £ RegiS* ,T*) 

(6) S) = bMh-j and T/ = V 7 e Unreg(5* , T*) 

We are now prepared to show that for generic parameters, (l**T i has exactly 2" fixed points. We will show 
that for an arbitrary set i? C [0 , « ) , there is exacdy one fixed point (5* , T * ) with Reg iS*,T*)=R. If R^Q), 
we have seen that {S*,T*) = (0,0). Under the assumption that R^%,S*: and TJ are determined for j G /? by 
(|5]l. Observe that if n is generic, then the T-values determined by (jSj are non-zero, so that Reg(5'*, T*) = R 
as required. Since the population values for j ^R are now determined by ^, we see that there is exactly 
one fixed point for each R C [0,«) as required. 

In general, not all of these non-vanishing fixed points (defined by (|5]l and (|6]l under the genericity as- 
sumptions) are biologically relevant, since a given pattern of regulated and unregulated stages can give 
Tj < for one or more values of j. We will call non-vanishing fixed points (i.e., different from the unin- 
fected equilibrium) infected equilibria. Infected equilibria with the property T* > V y = 0, ...,« — 1 (S* > 
V y = 0,...,« — 1 follows from the assumed positivity of parameters) are called biologically meaningful 
infected equilibria. 

4. Establishing infection: Stability of the uninfected equilibrium 

We start by deriving the expression for the basic reproductive number Rq already revealed above (see 
(O). Our first genericity assumption (O allows us to eliminate the case Ro= I. We will then show that the 
remaining dichotomy Rq < I 01 Rq > I determines whether the pathogen is able to establish infection. We 
give two proofs here. The first uses the interpretation of the parameters as describing the life-cycle of the 
pathogen. The second is based on the analysis of eigenvalues of the Jacobian at the uninfected fixed point. 
These two approaches correspond roughly to Rq and tq. For an excellent discussion of these, see 15|]. 

Given n and an arbitrary biologically meaningful point iS,T) = iSo,... ,Sn^i,TQ,... ,r„_i), we define 

Mj^ 'J^ 

oj+i+fj+i+Tj+i 

R = YlMj 

./=o 

Note that at a fixed point iS*,T*),Mj= Mj for j e Unreg(5*, T*). 
Proposition 1. 
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(1) The basic reproductive number of the system is Rq. 

(2) The reproductive number at an arbitrary biologically meaningful point [S, T) is R. 

(3) The overall reproductive number of the pathogen in the presence of any (non-negative) immune 
response is less than its basic reproductive number in the absence of an immune response. 

(4) If Rq < 1, the pathogen fails to establish infection and the uninfected fixed point (0,0) is a global 
attractor 

(5) At an infected fixed point, the reproductive number is exactly 1. Thus, at an infected fixed point we 
haveWj^p,,)Mj^l. 

Proof. To see that Rq is the basic reproductive constant of this system, consider the behavior of stage Sq when 
the naive host encounters one unit of this stage. The average Hfespan at this stage is l/(ao + /o)- During 
the course of that hfespan, it produces ^o/o/(ao+/o) units of stage Si. This, in turn produces {r\f\/{ai + 
/i))(ro/o/(flo+/o)) units of stage ^2 . Continuing in this way, we see that the original introduction of a unit 
of 5o resuhs in n"lo(o/;/(«i + /;)) = n%lMj=Ro units of 5o. 

To see that the basic reproductive number at an arbitrary point {S,T) is Wje\id,n)^i^ observe that in the 
presence of an immune response, the average Hfespan of stage Sj is l/{aj + fj + Tj). The result now follows 
as before. The third and fourth statements now follow immediately, while the final statement follows from 
observing that at a fixed point the pathogen exactly reproduces itself and therefore has basic reproductive 
number 1 . D 

We now examine Rq from a slightly different perspective. 

Proposition 2. IfRo < 1, the eigenvalues of the Jacobian matrix of system ^**| ) evaluated at {S,T) = (0,0) 
have negative real part and, consequently, the uninfected fixed point is asymptotically stable. On the other 
hand, if Rq > I, at least one eigenvalue of the Jacobian has positive real part and thus, the uninfected 
equilibrium is unstable. 

Proof. The characteristic polynomial P of J{S, T) is given by the following expression. 

P{X)^{~\rYi{{Sj-ai-X){b + X)-Sj{-{a, + fj) + b))-Y{{pi{b + X-Sj+i)) 

j=0 7=0 

where we use the notation Pj = rjfj and CT, = aj+fj + Tj. (See Theorem |5] of the Appendix for the details 
on how these expressions were derived.) Evaluating these at {S*,T*) = (0,0) gives 

P{X) = {-irib + Xr('fi{aj+fj+X)-'Yir,fX 

\j=Q k=0 ) 

Thus, the eigenvalues of /(0,0) are given by the «-fold root A = —b and the zeros of the polynomial/? (A) : = 

n— 1 ;;— 1 « — 1 

n (fli + fj+X)— n '"*:/*:• By our assumptions, all coefficients of R are positive except for H {'^j + fj) ~ 

;=0 k=0 ;=0 

n-1 

n '"aAi which is positive if and only if 

n nfk 

n{aj + fj) 
;=o 
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By Routh's criterion (yj]), if Rq> 1 , at least one root of R (and therefore of the characteristic polynomial P) 
has positive real part. Therefore, the equilibrium point {S,T) ~ (0,0) becomes unstable foiRo > 1. 

We now wish to see that if Rq < I, all roots of P have strictly negative real part implying that the equilib- 
rium point {S,T) = (0,0) is (locally) asymptotically stable. 

To see this, consider the polynomials 

n-l 






and 



E{X):=-llrjfj. 

j=o 

Thus, R = Q + E, and we can interpret R as a perturbation of Qhy E. By our assumptions, the real parts of 

all roots of Q are strictly negative, in other words, Q is strictly Hurwit20. Here the question arises as to how 

big the perturbation E can be, so that R remains strictly Hurwitz. By Corollary 4 in [|13il . this will occur i(j 



max 

(U>0 



E{ico) 



Qi'co) 



< 1 



where i is the imaginary unit. We have 

\E{io})/Q{ico}\-- 



-linfA/ilii^j+fj + i^) 

7=0 / \j=0 



'n-l 



\ I n nfj 1 / ( n(«; + fj + '■«) («7 + /; - '•«) 



fn-l 



'n-l 



Mino/.i /[Uii^j+fj)'- 



co^ 



It is now clear that this achieves its maximum, namely Rq, at co ~0. In particular, Rq < I implies that all 
eigenvalues of the Jacobian have negative real part and the uninfected fixed point is asymptotically stable. 

D 



5. Starvation: A partial order on the stages 

In this section, we introduce a relation j y k which we read as j starves k. This relation depends only on 
n. We will show that this relation limits which patterns of regulation can appear in a biologically meaningful 
fixed point, and will show that j y k is a strict partial order 

Definition 1. We will say stage j starves stage k and write j )~ k iffMji^ < 1 . 

The following proposition explains our use of the term starves. 

Proposition 3. Let the set of parameters it be generic. Furthermore, let {S*,T*) be a biologically meaning- 
ful infected equilibrium (in particular, Reg(5*, T*) ^ (d) and let j,k G [0,n) such that j ^ k. If j >- k , then 
j e Reg(5*, r*) implies k £ Unreg(5*, T*). 



A real polynomial is called strictly Hurwitz if the real part of each of its (complex) roots is strictly negative. 
We will see that such a maximum exists. 
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Proof. Suppose j G Reg(5'*, T*) and assume to the contrary that k G Reg(5*, T*). Let m G Nq be the number 
of regulated stages between j and k. That is, suppose [y + 1,A:) nReg(5*,r*) = {y'i,---,7m} and suppose 
these are ordered so that (cyclically) j — jo < j\ < 72 < •■• < jm < jm+i = k. Since {S*,T*) is biologically 
meaningful and ;ris generic, T* >0for7' GReg(5*,r*). By observation|3]of Section[3j wehaveM^^y^^j > 1 
for s ^0,... ,m. But this gives 

m 

contradicting the assumption that j >- k. D 

This result can be biologically interpreted as follows: If j >~ k, regulation at stage Sj starves stage S^ of 
sufficient population to support regulation. 

Proposition 4. Let % he such that Rq> I. Then y is a strict partial order 

Proof. We must show that >- is transitive, anti-reflexive and asymmetric. Suppose j )~ k and k )~ £. If 
k£ [j,£] then 

{Mj . . .Mk-x){Mk . . .Me_i) ^Mj.. .M^^i. 
Since each of the quantities on the left are less than one, so is their product as required. On the other hand, 
if £g [j,k], then 

(Mj . . .Mk_i)iMk . . .Me_i) ^Mj.. .Mj^iMj . ..Mt_, 

= RoMj...Mi^i < 1 

Since /?o > 1> 

Mj...Mi:^i<l, 

as required. 

To see that >- is asymmetric, note that R^ can be written as 

Ro^Mo.. .M„_i = {Mj . . .Mi_i ) {Mk . . .Mj^i) 

Thus Rq > I and j y k imply 1 < M^ . . .M,_i contradicting k>~ j. 

Finally, to see that >- is anti -reflexive, note that by definition Mjj = 1 . D 

Definition 2. We will say that Sj is unstarvable if j is >- -maximal. Otherwise Sj is starvable . We will refer 
to the set of indices of unstarvable respectively starvable stages as Str{7t) respectively Unstr{n). 

Remark 1. Note that if % is such that no two stages are comparable, then >- is empty and every stage is >- 
-maximal, so Unstr(7r) = [0,n). In particular, Unstr(7r) ^ V a:. 

6. Stability of infected equilibria 

In this section, we consider biologically meaningful infected fixed points. In particular, we can assume 
Rq > I- We will show that for generic n there is exactly one asymptotically stable biologically meaningful 
fixed point. We give two characterizations of this fixed point, one in terms of pathogen populations and one 
in terms of immune response. 

Definition 3. Let {S* ,T*) be an infected fixed point. We will say that the pathogen populations are moder- 
ated at (5*, r*) (or simply that {S*,T*) is moderated) ifS* < b for j G Unreg(5*, T*). We will say that the 
immune response is saturated at {S*,T*) (or simply that {S*,T*) is saturated) ifReg{S*,T*) = Unstr(;r). 

Remark 2. Given a fixed point (S*,T*), S* = bfor j G Reg(5*, T*). Thus at a moderated fixed point, at no 
stage is the pathogen population greater than the regulated population. 
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Proposition 5. Let the set of parameters % be generic. Furthermore, let {S*,T*) be a biologically meaning- 
ful infected fixed point. Then the pathogen populations are moderated at {S* ,T*) if and only if the immune 
response is saturated at (5*, T*). 

Proof. We start by proving that saturation implies moderation. Suppose (5'*,r*) is saturated, and suppose 
n is generic. If Unreg(5'*,r*) = 0, the claim holds vacuously. So suppose k G Unreg(5'*,r*). We then 
have A: S Str(7r). Letjbe ;^-maximal such that y )^ ^. lf[j+\,k) C Unreg(5'*,r*), we are done, for then 
SI = S*Mji, = bMjk < b. On the other hand if [; + 1 , yt) n Reg (5* , T* ) ^ 0, choose m e[j + l,k)n Reg (5* , T* ) 
so that [ot+1,^) cUnreg(5*,r*). Since M^^^ =Mj,„M,„i^ < 1 , we must either have M,,,,;. < 1 forcing 5^ <b, 
or we have M^,,, < 1. But the latter is impossible since m £ Reg(5'*, T*) and hence m G Unstr(7r). Hence 
{S*,T*) is moderated as required. 

Now we show the converse. Suppose {S*,T*) is moderated. We then have 

RegiS*,T*)^{j\S*^b} 

\JnKgiS*,T*) = {j\S*<b} 

We claim that Unreg(5*, T*) C Str(7r). If Unreg(5*, T*) = 0, this holds vacuously. Now assume Unreg(5*, T*) ^ 
0andletA:GUnreg(5*,r*). Choose j G Reg (5*, T*) such that [y + l,yt] cUnreg(5*,r*). ThenSl^S*Mjk = 
bMjif < b proving j y k as required. 

We claim that Reg(5*,r*) C Unstr(7r). Suppose to the contrary that k G Reg(5*,r*) nStr(7r). Then 
there is j G Unstr(;r) such that j >~ k. Since Unreg(5*, T*) C Str(7r), we must have j G Re.g{S*,T*). But if 
j G Reg(5*, r*) and j >- k, then, by Proposition^ k G Unreg(5*, T*) conti-adicting k G Reg(5*, T*). 

Since Unreg(5*,r*) C Str(7r) andReg(5*,r*) C Unstr(7r) and Str(7r), Unstr(7r) as well as Reg(5*,r*), 
Unreg(5*, r*) are partitions of [0,«), we have Unreg(5*, T*) = Str(7r) and Reg(5*, T*) = Unstr(7r) proving 
that {S*,T*) is saturated as required. D 

We now justify the use of the term "saturated". 

Theorem 1. Let the set of parameters % be generic. Furthermore, suppose {S* ,T*) is a biologically mean- 
ingful infected fixed point and assume {S* , T*) is not saturated. Then there is a j £ [0,n) so that T* = 
and for any neighborhood U of [S* ,T*), there is a biologically meaningful point x £ U so that Tj |^ > 0. In 
particular, {S*,T*) is unstable. 

Proof. Since (5'*,r*) is not saturated, by Proposition|5] it is not moderated. Furthermore, Unreg(5'*,r*) ^^0. 
Were Unreg(5*, T*) = 0, then the fixed point would be completely regulated (Reg(5'*, T*) = [0,n)) and such 
a fixed point is always saturated. To see this, assume that {S*,T*) is not saturated. Since Unreg(5'*, T*) = 0, 
this assumption is equivalent to Str(7r) ^ 0. Thus, there is a ^ G Str(7r) which is dominated by another 
stage j, in other words, j >- k. But j G Reg(5'*,r*) and therefore, by Proposition [3] k G Unreg(5'*,r*), a 
contradiction. 

Summarizing, there is j G Unreg(5'*, T*) so that S* > b. Let ctj be the unit vector in the Tj direction. 

Recall that ^{S*,T*) = {S* - b). Thus, for any sufficiently small 5 > 0, we have tjlp+Sej. > 0- This 
implies that arbitrarily close to {S*,T*), there are orbits which move away from {S*,T*). D 

Remark 3. One can think of (S* ,T*) as a state of a "micro-ecology" in which Tj plays the role of a species 
which is capable of invading. 

Theorem 2. Let (S*,T*) be a biologically meaningful infected fixed point. Assume that not all Tj are 
equal. If (5*, T*) is moderated then {S* , T*) is a local asymptotically stable equilibrium. In particular, the 
eigenvalues of the Jacobian matrix J{S* , T*) have strictly negative real part. 
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Proof. For brevity, we will use the notations M = Reg(5*,r*), '^ = Unreg(5*,r*), Pj = rjfj and (7/ = 
Qj + fj + T* (thus, Oj = cij + fj for all ; e Unreg(5*, T*)). The characteristic polynomial of J{S*,T*) is 
given by 

P{X) = ±li{{S*-aj-X)ib + X)-S*{-{aj+fj) + b))-fl{pj_,{b + X-S*)) 

.7=0 7=1 

By grouping the factors containing regulated resp. unregulated stages we obtain 

P{X) = ± n (A' + CT;A + T;b) n (A' + {aj +fj+b- S*)X + {aj + fj){b - S*)) 

-^'■^'n(p.-on(p.-i(^+^-^;)) 

je .«" je"?/ 

As in SectionlH we consider the polynomials 

(7) Q{X) : =±Y\{X^ + GjX + T*b)Y\{X^ + {aj + fj+b-S*)X + {aj + fj){b-S*j)) 

E{X) : =-Al-"^ln(p.--i)n(Pi-i(^+^-^}))- 
j&jf je'-?/ 

Thus, P = Q + E, and we can interpret P as a perturbation of Q by £. By our assumptions on the parameters, 
by the positivity of all T* for j G ^, and by the assumption S*; < b for all y G ^, Q is strictly Hurwitz 
(this is also correct in the case '^ = 0). Here the question arises as to how big the perturbation E can 
be, so that P remains strictly Hurwitz. By Corollary 4 in ||l3'|, the strict Hurwitz property is preserved if 
max(|£(/a))/e(/to)|) < 1.0 We have 



\E{ico)/Q{ico) 





('«)'-*i n (P;-i) n (pj-iiico+b-s*)) 




n (ii(oy + i(Jj(o 
je^^ 


+ T*b) n ({icof- + i{aj+fi + b-S*)co + 
' je"?/ ^ 


{aj+fj){b-S*)) 



7=0 je "?/ 



CO'' 



'S*)' 



n ^J(T*b-co^y + ajco'^ n \l({aj + fj){b^S*)-co^y + {aj + fj+b-S*y^co^ 



OJl' 



i"n(P;) n Jco'^ + {b^s*y 

;=0 je 9/ V 



n coJ{T*b-co^ /co^ + ajn J{{b-S*y- + (o^){(aj+fjY + (o^) 
jG # V V ■' / ^ je"?/ V ■' 

n iPj) 

;=o 



n \(T*b-(oA /co2+(t2 n 7(^7+77 



''"It is easy to verify that |£(0)/g(0)| = and Xmii^^^ \E{ia)/Q(i(o)\ = 0. Notice tliat Q(ia) is non-zero on (O,"). Since h((o) := 
\E{i(o)lQ(i(o)\ is a non-negative continuous function on [0, «>) (even differentiable on (0, oo)) and obviously not constant, it must take 
its maximum in tlie open interval (0,00). 
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Eachof the functions /i,( ft)) := (T*b- (O^j /(O^, j e^has a global minimum on (0,°°) at ft) = jTJb >0 

whereas the function g{(o) := {oj + fj)^ + (O^ has a global minimum on [0,oo) at ft) = 0. Therefore, we can 
conclude that for maX(B>o \E{i(o)/Q{ico)\, if existent, it must hold 

n-l 

nipj) 

max|£'(ift))/e(ift))| < ^ 






This is also coiTect in the case f/ = 0, in which \E{io})/Q{ico)\ = H {Pj) / I\ \ (T*b - (oA /ft)2 + aj 

7=0 je .«' V V ■' / -^ 

and Ylje"?/ i^i '^fj) ~ ^ • The strict inequality arises from the fact that not all T* are equal (and consequently 
not all \ T*b are equal) and thus, the functions hj{co) do not achieve their minima simultaneously. By 5. in 



Proposition [T] it follows maX(i)>o \E{ico)/Q{i(o)\ < 1 and the claim follows by Corollary 4 in 11 1311 . D 



Corollary 1. Let the set of parameters % be generic. Furthermore, let {S*,T*) be a biologically meaningful 
infected fixed point. Assume that not all Tj are equal. Then {S*,T*) is locally asymptotically stable if and 
only if{S*.,T*) is moderated. 

Proof. The claim follows immediately from Proposition |5] and Theorems [T] and |2l D 

The previous Corollary can be restated as follows. 

Theorem 3. Let the set of parameters K be generic. Furthermore, let [S* ,T*) be a biologically meaningful 
infected fixed point. Then [S* , T*) is locally asymptotically stable if and only if [S* , T*) is saturated. 

Proof. Under the assumption that n is generic, (5'*,r*) is saturated if and only if it is moderated. Note 
that for a generic parameter set K, the T* are not all equal as this occurs only on a subset of parameters of 
co-dimension n — 1 . D 

Remark 4. We have introduced a third requirement for % to be generic, namely that T* are not all equal. 
As we will see below, this is equivalent to the requirement that '";-i/;-i — {aj + fj) is not independent of j. 

Definition 4. Let {S*,T*) be a fixed point. IfReg{S*,T*) ~ [0,«), we will say that {S*,T*) is completely 
regulated. 

Remark 5. If {S*,T*) is completely regulated, it follows that [S* ,T*) is infected. It is also vacuously 
moderated. IfK is generic and [S* , T* ) biologically meaningful, it is also saturated (see the Proof of Theorem 
[7]). On the other hand, if K is generic and such that no two stages are comparable, then the completely 
regulated (not necessarily biologically meaningful) fixed point is saturated (cf RemarkU}. 

Corollary 2. Let {S*,T*) be the completely regulated biologically meaningful fixed point. Suppose that 
the Tj' are not all equal. Then (5*, T*) is locally asymptotically stable and the eigenvalues of the Jacobian 
matrix J{S*,T*) have strictly negative real part. 

Proof. Since {S*,T*) is completely regulated, {S*,T*) is moderated. D 



12 EDGAR DELGADO-ECKERT AND MICHAEL SHAPIRO 

Remark 6. Note that if we drop the assumption that not all the T* are equal, we are forced into the case 
T* = 1 for all i G [0, n)for some X G M.^ and consequently Unreg(5*, T*) ~ 0. In this case. 



( "-1 \ 



max|£'(;co)/<2('ft))| = max 

fi)>0 laX) 



n {Pi) 



«-l 



n-1 



nJirb-co'^r/co^ + sj 



n iPj) 



n-l 



VAoV^ ' ' J J -0 



Uiaj+fi + T*) 



and by Proposition\l] maXft,>o \E{ia))/Q{iCO) \ ~ I- Corollary 4 in 4731/ and the continuity of the eigenvalues 
on the entries of a matrix let us conclude that Re(A^) < Ofor all roots Xk of the characteristic polynomial P. 
In this very non-generic case, linearization of the right hand side of i\**\l does not yield a stability statement 
for (S*,T*). In the next subsection we will study this case in more detail. 

7. An embedded Lotka-Volterra system 

For certain highly non-generic sets of parameters, this system contains an embedded Lotka-Volterra 
predator-prey system. Recall that this system is given by 

x= {a ~y)x 
y={x-^)y 

where a,)3 G M^ are positive parameters and x and y are non-negative. (See, e.g., ij^.) We will use the 
notation A^ x Aj to denote points (^o , . . . , ^n _ i , Tq , . . . , r„ _ i ) where the values Sj and Tj are independent of 
j, that is, points satisfying ^o = ^i = • • • ~ Sn-\ and T(, = T\ = • ■ ■ — r„_i . 

Theorem 4. Suppose {S*,T*) is the completely regulated infected fixed point of (f**t. Then the following 
are equivalent: 

(1) Tj" is independent of j. 

(2) {S*,T*)eAsxAT. 

(3) The expression /"y-i//-! — {aj + fj) is positive and independent of j. 

(4) As X Aj- is invariant and the dynamics of the system restricted to Aj x Aj is the Lotka-Volterra 
system. 

(5) There is an embedded Lotka-Volterra system. 

(6) The Jacobian of the system has a pair of pure imaginary eigenvalues. 
Proof. 

m =^ |2]: We have T* ^ for < y < «, so S* = /? for < ; < n. 
|2] =^12: Recall that 

Since S* = b for < j <n and T* is independent of j , the result follows. 

13 =^ H]; We take Sj = x, Tj = y on As x At, O-1//-1 ^ («)+//) = <^ ^"^ b = ji.We then have 

Sj = n-ifj-iSj-i - {oj+fj + Tj)Sj 

= {a-y)x 

tj^{Sj-b)Tj 

^{x^P)y 

Since the right-hand sides are independent of j, we may replace the left-hand sides with the derivatives of x 
and y giving the Lotka-Volterra equations. 
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|5] immediately implies \5\ 

|6] now follows from the well-known dynamics of the Lotka- Vol terra equations. 

It remains to prove that|6] implies[7] This follows from Corollary|2] D 

Remark 7. The Lotka-Volterra dynamics take place on an embedded 2-plane. Corollary\2\says that the 
eigenvalues of the Jacobian have non-positive real parts. This implies that there is no exponential departure 
from the fixed point orthogonal to the Lotka-Volterra plane. We conjecture that the pure imaginary eigen- 
values of the Lotka-Volterra dynamics are the only ones with non-negative real part. If this conjecture is 
correct, the Lotka-Volterra plane is an attractor in a neighborhood of the fixed point. This conjecture is 
consistent with the results of numerical investigations. 

8. Self-establishing stages 

In this section we show how to modify Definition [T| to generalize Theorems [TIIJ] and [3] in the case where 
there exists a y £ [0,n) such that aj+fj < 0. 

In previous sections, we have assumed that aj+fj > for all j G [0,«). The assumption that a,- > says 
that stage j decays rather than proliferates. If aj < and aj+fj > 0, then stage j proliferates, but does so 
more slowly than it is lost to stage j+l. Accordingly, no stage is able to establish itself independently of the 
other stages. We now Uft this assumption and allow stages j such that aj+fj < 0. Such a stage proliferates 
faster than it differentiates and is thus able to establish itself independently of the remaining stages. We call 
such a stage self-establishing. 

Since fj and rj are positive for all j, a self-establishing stage also establishes infection at all other stages. 
The converse is also true. Suppose a naive is infected with a small amount of Sj > where j ^ SE(7r). 
Suppose also that k e SE(7r) with [j+l,k)n SE(7r) = 0. Given that Sj > for f = 0, (Q ensures that 5*: > 
for any sufficiently small t. S^ then proceeds to proliferate. 

As our fourth and last genericity assumption, we assume aj+fj 7^ V y £ [0,n). We take SE(7r) = {j \ 
aj + fj < 0}. In this section, we take it as a standing assumption that SE(;r) j^ 0. 

Proposition 6. Let nbe a set of parameters. Furthermore, let [S* ,T*) be a biologically meaningful infected 
fixed point. If j € SE(7r), then j e Re.g{S*,T*). 

Proof. If, on the contrary, T* = 0, it follows that Sj > 0, contradicting the assumption that {S*,T*) is a fixed 
point. D 

It should be clear that any self-establishing stage cannot be starved of sufficient population to support 
regulation. 

Definition 5. We say that j starves k and write j >- k ifMj^ < 1 and [7 + 1 , ^] fl SE(7r) =0. 

Note that this reduces to the previous definition when SE(7r) is empty. 

Proposition 7. Let it be such that Ro > 1. Then the relation >- is a strict partial order 

Proof. Let SE(7r) consist of ji < J2 < ■ ■ ■ < jm- These decompose [0,«) into Ui,...,Um+i where U( = 
Ijeje+i) for ^ — 1, ■ ■ ■ ,?« and U,„+i = [j„i,ji)- We order each of these segments according to the reverse of 
the order induced by the cyclic order of [0,«). That is, for£= I,..., m,Ue is ordered j /;> j ( + l> ... >j(+i — 
1, while f/m+i is ordered j„,> . . . >n — 1>0> . . . >ji — 1. Notice that if j >~ k then j and k lie in the same U( 
and j>k. It is easy to see that this makes >- a strict partial order. D 
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We define Str(7r), Unstr(7r), saturated and moderated as before. The statement of Proposition |3] as well 
as the equivalence of moderation and saturation (Proposition|5]l and their necessity for stability (Theorem[T) 
can be proved analogously, where the result of Proposition|6]plays an important role. 

In order to prove sufficiency in the presence of self-establishing stages, we need the following lemma. 

Lemma 1. Let the set of parameters % be generic. Furthermore, let (S*,T*) be a biologically meaningful 
infected fixed point. Then 

(1) S*>Qforje[Q,n). 

(2) Cj = aj + fj + T* > Ofor je[Q,n). 

Proof. Since {S*,T*) is infected and biologically meaningful, there is j so that S*; > 0. Hence, by ([T) 
ry/,-5* >0, whichinturnrequires -(fly+i+/,+i+r/_^i)5*^i < 0. This requires 5*^ j /O, and since (5*, T*) 
is biologically meaningful, S*^^ > 0. This, in turn, requires c7,+i = aj+\ +//+i + T*^^ > 0. Continuing in 
this way extends this to all j e [0,n). D 

The proof of Theorem |2] depends on equation (|7]), and we now revisit this in the case of a saturated 
fixed point with self-establishing stages. To see that 2(A) is Hurwitz, we must check that its coefficients 
are positive. For the product taken over regulated stages, we note that Oj and Tj'b are both positive. We 
turn to the product taken over unregulated stages. As we have seen, an unregulated stage cannot be self- 
establishing. In particular, at any unregulated stage, aj+fj is positive and we proceed as before. Theorem|3] 
and Corollary |2]now follow as before. 

9. Discussion 

In this paper, we have considered an ODE model i**\ of the interactions between a host and a pathogen 
which uses a cycle of antigenically distinct stages to establish and maintain infection. We were able to give 
a simple expression for Rq which determines whether the pathogen is able to establish infection. If Rq < 1, 
exposure to the pathogen, even in large quantities, will not succeed in infecting the host since the pathogen 
dies faster than it can replace itself. Indeed, the host will clear the pathogen even in the absence of an 
immune response. 

If Rq > 1, then the pathogen is capable of infecting the host, and in this model, an immune response is 
necessary to control the infectiono In this case, the absence of an immune response leads to unbounded 
expansion of the pathogen population and presumably, the host's death. 

The question then arises, when an immunocompetent host is infected with an infectively competent 
pathogen (i.e., Ro> 1), can they arrive at a stable equilibrium? If so, what are the balancing populations for 
pathogen population and host response? In terms of the model, this question translates to asking which of the 
model's equilibria are stable and under what conditions. A key ingredient towards answering this question 
was the ability to distinguish the different stages of the pathogen in terms of the net yield with which they 
are generated. We were able to formalize this distinction using a binary relation h >~ k, which turns out to be 
a partial order if Ro> 1 (see SectionlSj. 

Assuming the pathogen can establish infection the concepts of starvable and unstarvable stages restrict 
the possible patterns of regulation. The starvation relation j y kis based on comparison of the rate at which 
stage j produces stage k with the rate at which stage k is lost to death and differentiation. This determines 
the ability of stage j to replenish the population at stage k. When this downstream yield is smaller than the 
regulated population b, it is insufficient to support regulation. In this way, if j >- k and j is regulated, then 
regulation of stage k vanishes. In this sense, the immune regulation of unstarvable stages is sufficient to 



In models which include the supply of uninfected tissue, the rate of supply provides a limiting factor for infection in the absence 
of immune response. 
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immunologically control the starvable stages. At steady state, immune regulation is only required against 
those stages that are produced with relatively higher yield and the equilibrium becomes moderated (see 
Section |6]l. Accordingly, there is a unique biologically meaningful infected equilibrium, namely, the one at 
which all unstarvable stages are regulated and all starvable stages are unregulated. 

If the parameters of the system allow for the existence of a biologically meaningful completely regulated 
fixed point, the model predicts that this will be the (only) stable equilibrium under the premise that the 
pathogen can establish infection (Corollary |2l). According to equations (|5]l, the completely regulated fixed 
point satisfies S* ^ b and T* = rj^ifj^i — (oj+fj), j = 0, ...,n — 1. The inequalities T* > imply that the 
amplification factors rj are big enough for every stage to be replenished at a higher rate than the sum Oj + fj 
of its own decay and transformation rate. Our model predicts that a pathogen which produces every stage 
with such an effectivity can only be immunologically controlled with a fine tuned immune response against 
every stage the pathogen cycles through. 

Having clarified the properties of equilibria, it is pertinent to keep in mind that these are local properties 
of the model. We were unable to perform a global analysis of the dynamics (see, for instance, [19]). That is, 
we cannot give an account of trajectories along which the system might approach such a steady state config- 
uration, the basins of attraction of attracting fixed points, nor can we eliminate the possibility that the system 
exhibits limiting cycles or other attractors. Our naive expectation is that the unique biologically meaningful 
infected equilibrium acts as a global attractor in the case where the host possesses an immunocompetent 
response to every stage. 

Our methods for establishing the uniqueness and existence of the unique stable fixed point have very dis- 
tinct flavors. The arguments for uniqueness have clear biological interpretations. In particular, the starvation 
relation, j >~ k of Section |5] considers the maximal rate at which stage j can produce stage k with the rate 
at which stage k is lost to stage A: + 1 and in turn compares these to the minimum population necessary to 
support immune regulation. Likewise, Theorem[T]describes the ability of an immune response to invade and 
therefore alter a given equilibrium. 

Our methods for establishing existence have a very different character. Here we have turned to the eigen- 
values of the system's Jacobian matrix and these have lead us to the roots of its characteristic polynomial. 
These sorts of computations can offer many challenges which have been the subject of intense study for the 
last five decades within the control and systems theory community. Many of the resulting techniques focus 
on the sorts of perturbations that can be applied to a strictly Hurwitz polynomial without losing the Hurwitz 
property. We predict that as biological models continue to grow in complexity, mathematical biologists will 
increasingly turn to these analytic tools. For entry into this topic we recommend the website |14]. 

Having stated the general predictions that can be derived from our model, we would like to elaborate on 
the model's features, properties and predictions within the biological framework of chronic EBV infection. 
It was indeed human EBV pathology what inspired the basic structure of our model. To this end, we provide 
a brief overview of what is known about EBV infection in humans. 



EBV is a highly successful pathogen which infects over 90% of the adult human population II22I1 . It is 
horizontally transmitted via saliva \W\ and is tropic to the epithelium of the oropharynx and B-cells. It infects 
naive B-cells and causes them to become activated B-blasts. These enter the germinal centers of Waldeyer's 
ring where they are ultimately transformed into latently infected quiescent memory B-cells. These circulate 
in the peripheral blood where they rarely express viral proteins and are therefore invisible to the immune 
system ||9[]. Upon return to Waldeyer's ring, some of these are triggered to start producing virus il2ll . These 



lytically infected B-cells ultimately burst, producing free virus which may be transmitted horizontally or 



may reinfect additional naive B-cells, thus completing the cycle 112411 . The epithelium of the oropharynx 
appears to play a role in amplifying virus for transmission \^. The host mounts T-cell responses against the 
blast, germinal center and lytic stages, and an antibody response against the free virus. The normal course of 
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acute infectious mononucleosis involves an acute phase during which latently infected memory B-cells can 
become as much as 50% of the memory B-cells in the peripheral circulation jlOtl . Immune response is also 
high during the acute phase. This resolves to a chronic phase of persistent low-level infection JlOll . Illlll . 

This suggests that a mathematical model of EB V infection might include the following infected popula- 
tions: activated B-blasts, infected germinal center B-cells, quiescent infected memory B-cells and lytically 
infected B-cells. These latter can be subdivided into stages exhibiting distinct antigenic profiles, namely 
immediate early, early and late antigens. Due to its short life span one might choose to omit free virus from 
such a model. In short, it is not at all obvious what would be the "right" number of stages to include in an 
EBV model, and this has motivated our choice to study a model which is entirely general as to the number 
of stages. 

We have chosen to omit the supply of uninfected tissue from our model for the sake of symmetry. We 
would argue that this does not affect the validity of the model unless a significant fraction of the uninfected 
tissue is lost to infection. To see this, consider a modification of the model to take this into account. Let us 
use S„ to model uninfected tissue, let us take Sn-i to model free viruo and let the rate of new infection be 
proportional to their product. We then have 

Sn-l = r„-2fn-2Sn-2 - (««-! +i85„ + T„^l)S„-l 

S„ = X- {a„ + r„j35„_i )5„ 

Sq = r„^S„S„-\ - {ao + fo + To)So 

(Compare equation 3.1 of 11611 .) If Sn is relatively constant, then S„ can be eliminated from the model and 
we can take /„_i = /35„ and r„_i/„_i = r„j35„. 

While acute phase EBV infection shows infected cells occupying a large portion of peripheral memory, 
it seems unlikely that it consumes a large portion of the uninfected naive B-cells supplied to Waldeyer's 
ringlM], 12111 . Accordingly, it seems unlikely that this will be an obstruction to application of this model. 

Our primary purpose in developing this model is to cast light on the long-term persistence of EBV infec- 
tion. Here we can surely treat the supply of uninfected tissue as nearly constant. 

Biologists will be interested in solutions to Q rather than (l**T i. The follow-on factors are unchanged. 
The analogs of equations ( |5]l and ^ are 




S* = - and r; = - ^r,_i/,_iM,^,_i - (a, + /,) V ; G Reg{S*,T*) 



S*j = —Mh.j and T/ = V y £ Unreg(5*, T*) 

Chj 

The starvation relation is now j )~ k if Mj/^ < -^. 

People show different patterns of infection and immune response to EBV. Under 10% of the adult pop- 



ulation is EBV negative. Given that over 90% of the adult population is EBV positive 112211 . 112011 . it is clear 
that most of the EBV negative population has had repeated exposure to EBV. Since these people show no 
EBV-specific antibodies |611, it seems clear that they are negative due to the virus's inability to establish 
infection rather than the immune system's ability to clear it. This shows that Rq is a host-pathogen property, 
not simply a property of the pathogen. 

Chronically infected people show very stable levels of infected memory B-cells in the peripheral blood 



I n|] . This is an indication that we are dealing with a stable fixed point of the host-pathogen system. However 



the level of infected memory B-cells varies from person to person Illlll as does the pattern of immune 



In this case 7^,_i is liumoral response. 
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response. To the best of our knowledge, there is always a humoral response to the free virus. Beyond that, 
we have seen the individuals who exhibit T-cell responses to both the blast and germinal center stages and 
to the lytic stages and individuals who only exhibit detectable T-cell response to the lytic stages |2]. The 
latter patterns could arise from a gap in the available T-cell repertory. We think it is at least as likely that the 
starvation relationship varies due to population variation in K. Empirical determination of these parameters 
is work in progress. 

It is instructive to consider the variation levels of latently infected memory B-cells. To the best of our 
knowledge, this compartment is never regulated llOll . Due to the very low levels of viral protein expression, 
this stage has a very low overall antigenicity as compared to any other stage. Taking S'latent to be latently 
infected memory B-cells, we see that a low value for ciatent compared to other stages implies starvability 
since we then have M.iatent < -^^■ 

■' Hatent 

The expression S* = ^, shows that the size of regulated populations depends on factors governing im- 
mune response, not on the flow of cells into that population. In particular, we can expect the size of the 
lytically infected B-cell population to show little or no correlation with the size of the latently infected 
memory B-cell population, while it might well be a correlate with immune factors such as HLA types. 

The question arises whether to expect a correlation between Tj* j^, and S'jtjtenf Here the picture is not clear 
as the model's predictions depend on the reasons for the variation in 5'*[jj(g|jj. If high values of Sl^^^^^ are due to 
a high rate flow into this compartment from the infected germinal center population, then we would expect 
a positive correlation between Tj* j^, and 5']*jjtent- ^f' ^'^ ^^^ other hand, high values of 5'(*3jgjjj are due to low rate 
of flow from the latently infected memory B-cell compartment into the lytically infected population, then 
we would expect a negative correlation between Tj* j^ and ^jt^j^j^j. 

On a different vein, from the point of view of the treatment of chronic infections, the characteristics of 
the immune response revealed by our model, in particular, the possibility of unregulated stages, give rise to 
speculation on the possible targets of such a treatment. Which stages would be most suitable drug targets, 
regulated or unregulated ones? 

Finally, we would like say a few words about some mathematical questions raised by this model. Our 
stability results are all local. We would conjecture that the unique asymptotically stable biologically mean- 
ingful infected fixed point is a global attractor for the positive orthant. This would follow from the existence 
of an appropriate Lyapunov function. However, we have not yet managed to find one. 

We also conjecture that in the Lotka-Volterra case, the 2-plane that carries the Lotka-Volterra dynamics 
is an attractor, at least in a neighborhood of the fixed point. Our numerical investigations to date support this 
conjecture. 

There are a number of generalizations of this model which deserve investigation. These include the supply 
of uninfected tissue, and immune cross-reactivity in which an immune response might target multiple stages. 

One interesting generalization concerns the situation in which each stage may be regulated by an immune 
response to more than one epitope and in which each of these epitopes occurs in multiple variants. This 
question is addressed for single-stage pathogens in JITT . Indeed Nowak suggested |15il that we examine this 
question for multi-stage pathogens and we have preliminary results in the case where only one stage exhibits 
this phenotype. 
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Appendix 
Derivation of an expression for tlie cliaracteristic polynomial of 7(5, J). 

Theorem 5. The characteristic polynomial of the Jacobian J(S, T) of system ( 1**1 ) is given by 



P{^) = i-irYliiSj-aj-X)ib + X)-Sji-{aj + fj) + b))-"Yl{Pjib + X-Sj+,)). 



Proof. The partial derivatives of the system are given by 



^ = < -fl/t - fk - Tk if j = k 
[ otherwise 

dTj |o otherwise 

dG^^hk if j = k 
dSj 1 otherwise 

dG^^iSk-b if j = k 
dTj 1 otherwise 



Thus, we have 



JiS,T)-(^C I 
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where A,B,C and D are the following nxn matrices 

/-C70 ■ 



A = 



B 



C 



D 



Pa -(y\ 



V ... 

/-5o 


V ... 

(To . . 

■■ 
Vo .. 

(Sa-b 




Pn-\ \ 





P„-2 -On-\j 

... \ 


-5„_i/ 




Tn-xj 



\ 




V ••• Sn-i-bj 

using the notation Pj = Vjfj and (7/ = a^ +// + 7}. It is well known that a similarity transform does not 
modify the eigenvalues of a matrix. The similarity transform we use consists of the following operations; 
The first n columns of 7(5, T) are transformed to 



The resulting matrix 

undergoes the following row transformation 

(C" D') 
The resulting matrix 



A' B 
C D 



(C D) + {A' B) 



'A' B 

then undergoes the following column permutations 

(B A' 

\D' C" 

followed by row permutations 

^D' C" 

B A' 
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The similar matrix obtained satisfies 

/So - do 



A'- 



PO Si - Gi 





Pn-l \ 







C": 



V ••• p„_2 S„-i-a„-iJ 

/To-Go + b ■■■ p„_i 

po Si-ai+h '■■ 

■•. ■. ■. : 



D' 



V 
f-b 

■. 





Pn-2 T„^l-(J„-l+bJ 







\ • • • -bj 
Now we look at the eigenvalue equations of the transformed matrix 



B A' 



D' C" 



B A' 



A -XI 



where / is the 2n x 2n unity matrix. If we assume X ^ —b, the matrix 

/-{b + X) ••• \ 

^'= ° '■ ' 

V •••0 -{b + X)J 

is invertible and we can use a well known formula for the determinant of block-partitioned matrices 



e(A):=det 



B A' 



det(D')det(A'-BD' C") 



[So — Gq — X —So 



{To-ao+b) 



det 



b+l 



„ _ Sop,i-[ 
Pn-\ b+X 






Pn-2-^-^^^ Sn-X-On-l-X-S,,- 



-1 b+X 
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Laplace's theorem applied to the last column of the matrix yields for A G C\{— /?} 

QiX) = i-ir(YiiiSj-aj-X){b + X)-Sji-{aj+fj) + b)) + {-ir+''Y\{pj{b + X-Sj^,))] 
\j=o j=o ) 

Summarizing, we have obtained an expression 2(A) that is equal to the characteristic polynomial P(A) 
of J{S,T) for all A G C\{— Z?}. If we drop the requirement A G C\{— /?}, Q is a polynomial on C and 
thus a holomorphic function on C As a consequence, we have two holomorphic functions Q • C — > C and 
P : C ^- C that are identical on the connected open subset C\{— ^}. By the well-known identity theorem for 
holomorphic functions, Q = P must hold. D 
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